Quantum limits in image processing 



V. Delaubert 1 ' 2 , N. Treps 1 , C. Fabre 1 , H.A. Bachor 2 , and P. Refregier 3 
1 Laboratoire Kastler Brossel, Universite Pierre et Marie Curie-Paris6, 
CC74, 4 Place Jussieu, 75252 Paris cedex 05, France 
2 ARC Centre of Excellence for Quantum-Atom Optics, 
Australian National University, Canberra, ACT 0200, Australia and 
3 Physics and Image processing group, Fresnel Institute, UMR 6133, 
Domaine Universitaire de Saint-Jerome, 13397 Marseille cedex 20, France 
(Dated: February 1, 2008) 

We determine the bound to the maximum achievable sensitivity in the estimation of a scalar 
parameter from the information contained in an optical image in the presence of quantum noise. 
This limit, based on the Cramer-Rao bound, is valid for any image processing protocol. It is 
calculated both in the case of a shot noise limited image and of a non-classical illumination. We 
also give practical experimental implementations allowing us to reach this absolute limit. 
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Introduction. - Many of the most sensitive techniques 
for the measurement of a physical parameter - which we 
will call p in the following and assume to be a scalar - 
are optical. In some cases, the total intensity or am- 
plitude of the light beam varies with p and conveys the 
information. It is then well known [HE El El that there 
exists a standard quantum limit in the sensitivity of the 
measurement of p when the light beam is in a coherent 
state, and that it is possible to go beyond this limit us- 
ing sub-Poissonian or squeezed light. In other cases, the 
parameter p of interest only modifies the distribution of 
light in the transverse plane and not its total intensity. 
The present paper deals with this latter situation. For 
example, the parameter p modifies the position or direc- 
tion of a light beam. This configuration has been studied 
at the quantum level, both theoretically and experimen- 
tally using split detectors [EIS0 or homodyne detection 
0, • But in many instances the parameter p affects in 
a complicated way the field distribution in the detection 
plane (that we will call here the image). For example a 
fluorescent nano-object imbedded in a biological environ- 
ment modifies the image recorded through a microscope 
in a complex way because of diffraction. Nevertheless 
its position can be determined from the information con- 
tained in the image with a sensitivity which can be much 
better than the wavelength 10] . In order to extract the 
parameter value in such experiments, one needs to use 
detector arrays or CCD cameras and to combine in an 
appropriate way the information recorded on the differ- 
ent pixels. 

When all the sources of technical noise have been re- 
moved in the apparatus, quantum fluctuations still affect 
the optical measurement and limit its sensitivity, in a 
way that can be readily calculated for each specific mea- 
surement protocol. The purpose of this paper is much 
broader. It is to answer the following question: what is 
the lowest limit imposed by quantum noise to the accu- 
racy of the determination of p, independently of the in- 
formation processing protocol used for the extraction of 
information? As we will see, this optimum limit depends 



only on the statistics of the fluctuations of the incom- 
ing light. We use an approach based on the Cramer-Rao 
Bound. This tool, widely used in the signal processing 
community |ll| , has already been applied to different do- 
mains, such as gravitational wave detection |T^ or dia- 
magnetism |13| . 

Notations and assumptions. - The parameter p is mea- 
sured relative to an a priori value chosen for simplicity to 
be 0. Because of the quantum fluctuations in the optical 
measurements, there will be an uncertainty on its esti- 
mation. An evaluation of this uncertainty thus provides 
the precision on the determination of the parameter p 
around a zero value. 

The mean value of the local complex electric field op- 
erator in the image plane, normalized to a number of 
photons, will be written for a given value of the param- 
eter as 

E{r,p) = 2VNu {r,p), (1) 

where N refers to the total number of photons detected in 
the mean field during the integration time of the detector. 
N is assumed to be independent of p, as stated previously. 
Uo(r,p) is the p-dependent transverse distribution of the 
mean field, complex in the general case, and normalized 
to 1 (its square modulus integrated over the transverse 
plane equals 1). The local mean photon number detected 
during the same time interval is 

n(r, P ) = ^^ P) =N\Mr,P)\ 2 - (2) 

Intensity measurements. - We first assume that p is 
determined by processing the information contained in 
the measurement of the local intensity, i.e. local number 
of photons. The best possible local intensity detection 
device would consist in a set of indexed pixels paving 
the entire transverse plane, in the limit when their spa- 
tial extension approaches 0. Let n = [rii, nk, ■■■] be 
one measurement of the photon distribution with such a 
hypothetical detector, where corresponds to the num- 
ber of photons detected on pixel k. Because of the noise 
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present in the light, the sample n differs from its statis- 
tical mean value n(p) = [n%(p), ...,n^(p), ...]. Let L(n\p) 
be the likelihood of its observation. Note that n cor- 
responds to a single measurement and hence does not 
explicitly depend on p, contrarily to the average on all 
the possible realizations. 

The achievable precision on the estimation of p is lim- 
ited by the Cramer-Rao Bound (CRB). More precisely, 
the variance of any unbiased estimator of p is necessarily 
greater than cr^ in = I/If, where the Fisher information 
If is given, when the actual value of p is 0, by [TT| : 



dp- 



■l(n\p) 



L(n|0)dn, 



(3) 



where we have introduced the log- likelihood l(n\p) = 
lnL(n|p). The integration spans continuously over all 
possible photon distributions that can be detected when 
the parameter value is p — 0. This information is thus 
highly dependent on the noise statistics. 

Let us first assume that the illumination is coherent, 
and therefore that the local intensity noise is Poissonian. 
The probability of measuring photons on pixel k, when 
the parameter equals p is given by 
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Restricting our analysis to spatially uncorrelated beams, 
the likelihood L(n\p) simply corresponds to the product 
of all local probabilities given in Eq. 0] Then, taking the 
limit of infinitely small pixels and using Eq. |3 one can 
show that the Fisher information equals 



jPoisson 
1 F 



n'(r,0) 2 
n(r, 0) 



n"(f,0) 



(5) 



where the ' denotes a derivative relative to p. Using Eq.^ 
one finally finds that 
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where a is a global positive parameter characterizing the 
variation of the image intensity with p, defined by 
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The smallest value of p that can be distinguished from 
the shot noise - i.e. corresponding to a signal to noise 
ratio (SNR) equal to 1 -, whatever the algorithm used to 
determine it from the local intensity measurements, pro- 
vided that it gives an unbiased estimation of p, is finally 
greater than a/2^/N. This value sets the standard quan- 
tum noise limit for intensity measurements of p, imposed 
by the random time arrival of photons on the detector. It 
is inversely proportional to the square root of the num- 
ber of photons, as expected, and is related, through the 



dependence with a, to the modification of the mean in- 
tensity profile with p. 

We now consider a non-classical illumination, still with 
identical mean intensity, but with local sub-Poissonian 
quantum fluctuations described by a noise variance 
up < 1 (assumed to be the same over the entire trans- 
verse plane). One can show that the CRB leads to 
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As we have already noticed, the limit given by Eq. |H1 
is valid for any measurement strategy. Nevertheless, a 
practical way enabling us to reach such an absolute limit 
remains to be found. This is what is presented in the 
next paragraph. 

Let us assume that an image processor calculates a 
given linear combination of the local intensity values 
recorded by the pixels of an array detector, as represented 
in Fig. H Assuming that the pixels are small compared 
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FIG. 1: Array detector as an optimal intensity detection, g 
refers to the pixel gain distribution. S(p) is the p-dependent 
processed signal. 



to the characteristic variation length of the image, the 
mean value S(p) of the computed signal can be written 
as an integral over the transverse plane as follows 




S (P)= / g(r)n{r,p)d 2 r, 



(9) 



where g(r) is the local gain on the pixel localized at posi- 
tion r, which can be positive or negative. Assuming that 
p is small, \uo(r, p)\ can be expanded at first order into 



\uo(r,p)\ = \u o {f,0)\ + -u/(r), 



(10) 



where a has been defined in Eq. \7\ and ui is a transverse 
function normalized to 1, already introduced in a partic- 
ular case in reference 0. The gains are chosen such as 
S(p = 0) = (difference measurement), so that S(p) is 
given, at first order, by 



S(P) 



2Np 



g(r)\u (f,0)\u 1 {r)d 2 r, 



(11) 



For a coherent illumination, because the local quantum 
fluctuations are uncorrelated, the noise variance A5* 2 on 
on S(p) i s eq ual to the shot noise on each pixel weighted 
by g 2 {r) [l4[- Moreover, as p is small, the noise is inde- 
pendent of p at the first order, and we get 



AS 2 = N 



! (r)\u Q (r,0)\'d'r 



(12) 
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It is then possible to optimize the gain factor g(f) 
in order to get the highest possible SNR defined by 
SNR= S(p) 2 /AS 2 . Using Cauchy-Schwartz inequality, 
one can show that the highest SNR value for the present 
measurement strategy is given by SNR= ANp 2 /a 2 , and 
is obtained for an optimal value of the gain distribution 
given by 



9opt{r) = (3- 



ui(r) 



(13) 



Mr, 0)|' 

where /3 is an arbitrary constant. The minimum measur- 
able value of p - corresponding to comparable signal and 
noise, i.e. SNR= 1 - is given by p m in — a/2yN, which 
is precisely equal to the CRB for classical illumination. 
The present measurement strategy is therefore optimal 
as it allows to reach the CRB for small values of the pa- 
rameter p, with the certainty that no other measurement 
strategy can do better. Note that we have not proven 
that it is the unique way to reach the CRB. 

We can even extend this result to the use of non clas- 
sical light. Indeed, using a bimodal field composed of a 
bright mode uq in a coherent state carrying the mean 
field, and a squeezed vacuum mode in the mode uj, the 
detected noise power is modified into 



AS* 2 = Nap, 



(14) 



when the noise variance on the amplitude quadrature of 
the mode ui is given by ap . Note that the use of a lo- 
cally squeezed beam would not provide any improvement 
as ui is the only mode contributing to the measurement 
noise, referred to as the noise-mode of detection [si Il4j|. 
The minimum measurable p value is in this case 
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(15) 



We have thus found a way to reach the bound, i.e. the 
minimum accessible p- value given in Eq|S] Moreover, our 
scheme requires minimum quantum resources, namely a 
bimodal field with squeezing in only one mode. 

Field measurements. - We now assume that the infor- 
mation about p is extracted from the knowledge of the 
local complex field, i.e. local amplitude and phase, that 
can be obtained by interferometric techniques. Similarly 
to the previous section, the best possible detection would 
here access the local complex field on k-indexed areas 
paving the entire transverse plane, in the limit when their 
spatial extension approaches 0. Let E = [E±, E k , •••] 
be a single measurement of the field distribution, hence 
independent of p, where Ek corresponds to the complex 
field detected on area k. Again, because of the noise 
present in the light, the sample E differs from its statisti- 
cal mean value E(p) = \Ei(p), ...,E k (p), ...]. Let L(E\p) 
be the likelihood of its observation. 

A bound to the maximal precision on p can again be 
calculated from the CRB, which, for a field measurement, 
is the inverse of the following Fisher information 



If = - 



L(E|0)rfE, 



(16) 



We assume that the local field fluctuations in the trans- 
verse plane can be described by a Gaussian probability 
density function independent of the mean field. More- 
over, we consider a classical or non classical illumination 
whose amplitude and phase quadrature fluctuations are 
described by the noise variances ap and <7q , respectively. 
These factors neither depend on r nor on p, as we assume 
the fluctuations to be homogeneous and independent of 
the parameter p. 

The probability to measure a field given by E k = P k + 
iQk on area k, where P k and Q k correspond to the local 
field quadratures, is, for a parameter value p 
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where P k (p) and Q k (p) are the local quadratures sta- 
tistical averages, satisfying: E k (p) = P k (p) + iQ k (p). 
Without loss of generality, we define the orientation of 
the Fresnel diagram relative to the phase of the mean 
field, i.e. P(f,p) = 2^/Nuo(r,p) and Q{r,p) — 0, taking 
the limit of infinitely small detection areas. Assuming 
no spatial correlations in the field fluctuations, one can 
show that the Fisher information is given by 



rGauss 
1 F 





\dP(r,p)l 


\ o 2 p 


dp 



d 2 r. 



(18) 
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We now introduce b, a second global positive parameter, 
characterizing the variation of the image field with p 



1 



d_ 

dp 



uo{r,p) 



d 2 r 



(19) 
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Using Eaniand ll9l the Fisher information simplifies into 
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The smallest value of p that can be distinguished from 
the quantum noise using a field detection of the optical 
beam is finally greater than the CRB: 
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(21) 



Again, we can propose an experimental scheme en- 
abling to reach this limit in the case of a small parameter 
p, and for which the mean value of the complex electric 
field can be written at first order 



E(r,p) = 2VN 



u (f,0) + -u E (r) 



(22) 



where b has been introduced in Ea ll9l and where ue(t) 
is a transverse function normalized to 1. 

Let us consider a balanced homodyne detection, as rep- 
resented in Fig. [21 with a local oscillator (LO) chosen to 
be defined by the following complex field operator 



p=0 



E LO {r) = 2^N^>u E {r)e 



,i0LO 



(23) 
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FIG. 2: Balanced homodyne detection as an optimal field 
detection. n-(p) is the p-dependent signal. (f>LO is the local 
oscillator phase. BS : 50/50 beamsplitter. 



Moreover, when the component of the image selected by 
the LO is in a non classical state , i.e. allowing a squeezed 
vacuum ue mode with a noise variance cxp on the ampli- 
tude quadrature within the incoming beam, we get 



bap 



Pmin — 



2VN 



(27) 



where N^o corresponds to the number of photons de- 
tected in the entire LO beam during the integration 
time. The LO is much more intense than the image, 
i.e. Nlo ^S> N. 8lo is the LO phase. The mean inten- 
sity difference n_ between the photocurrents of the two 
detectors is given in terms of number of photons by 

n_(p) = i J [E* LO (r)E(r,p) + E*(r,p)E LO (r)] d 2 r(2A) 

Note the similarity with Eq|5| as incident field amplitude 
and LO field play here identical roles of incident inten- 
sity and electronic gain, respectively. Though detectors 
do not resolve the spatial distribution of the beams and 
no processing of the spatial information is made, the bal- 
anced homodyne technique directly provides an " analog" 
computation of the quantity of interest. 

When the LO phase is tuned to the maximum of the 
p-dependent term, the homodyne signal becomes 

n_ = 2^NN LO j- (25) 
o 

For coherent illumination, the noise power on the homo- 
dyne signal corresponds to N^o, i-e. to the shot noise 
of the LO. The SNR of the homodyne measurement is 
thus given by SNR= 4Np 2 /b 2 . The minimum measur- 
able value of p - corresponding to a SNR of 1 - with 
homodyne detection is given by 

b . . 

p mm = — =. (26) 
2-//V 



This result corresponds exactly to the CRB calculated 
for amplitude measurements of p. The homodyne detec- 
tion scheme with the appropriate LO shape and phase is 
therefore an optimal field detection of p. Again, it uses 
minimal resources as only one source of squeezed light is 
needed to reach the non classical CRB. 

Comparison and conclusion. - We have presented two 
efficient - i.e. reaching the associated CRB - signal pro- 
cessing techniques for the extraction of information con- 
tained in an image. We can show that the CRB for field 
measurements, in which all quadratures can be accessed, 
is smaller than the one for intensity measurements, i.e. 
a > b. Yet, both schemes are useful: the intensity scheme 
is interesting since it is not restricted to monochromatic 
light, whereas the amplitude scheme is useful since it does 
not require pixellized detectors. 

Let us finally note that this work provides limits which 
are valid for a shot noise limited light of any shape. How- 
ever it is not valid so far for any kind of non-classical 
light, as we have restricted our analysis to homogeneous 
squeezed states. In a forthcoming publication, we will 
study the situation of quantum fields with quantum spa- 
tial correlations in the transverse plane and will investi- 
gate for the corresponding modifications of the CRB. 
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